Testing the Elliott- Yafet spin-relaxation mechanism in KCg; a model system of biased graphene 
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Temperature dependent electron spin resonance (ESR) measurements are reported on stage 1 potassium doped 
graphite, a model system of biased graphene. The ESR linewidth is nearly isotropic and although the g-factor 
has a sizeable anisotropy, its majority is shown to arise due to macroscopic magnetization. Albeit the homo- 
geneous ESR linewidth shows an unusual, non-linear temperature dependence, it appears to be proportional 
to the resistivity which is a quadratic function of the temperature. These observations suggests the validity of 
the Elliott- Yafet relaxation mechanism in KCg and allows to place KCg on the empirical Beuneu-Monod plot 
among ordinary elemental metals. 

PACS numbers: 76.30.Pk, 71.70.Ej, 75.76.+J 



I. INTRODUCTION 

Information storage and processing using electron spins, 
commonly referred to as spintronicsi, is an actively studied 
field. Spintronics utilizes the prolonged conservation of the 
spin quantum number, as the spin-relaxation time (Ti) usually 
dominates over the momentum relaxation time (r) by several 
orders of magnitude. In metals with inversion symmetry, Ti 
is described by the so-called Elliott- Yafet (EY)=i^ theory. The 
EY theory describes that the otherwise pure spin-up and down 
states of the conduction band are admixed due to spin-orbit 
coupling (SOC). The strength of the mixing is given by ■^, 
where L is the matrix element of the SOC for the conduction 
and a near lying band with an energy separation of A. The 
conduction spin states being admixed, spin-flip transitions are 
possible with a low probability whenever momentum scatter- 
ing occurs. Through first-order time dependent perturbation 
theory, Elliott connected the g-factor shift, the spin and the 
momentum relaxation times as: 
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where go ~ 2.0023 is the free electron g-factor, q;i,2 are con- 
stants around unity and are determined by the band structure. 
The former of these equations is known as the Elliott-relation. 
The EY theorji^ explained spin-relaxation for most monova- 
lent elemental metals^-- and later studies showed its validity 
in one-dimensional'' and polyvalent^'*^ metals and its gener- 
alization explained the spin-relaxation in metals with strong 
correlations-'". 

The recent discovery of graphen&i^ directed the attention of 
spintronics research towards carbon nanostructures. The weak 
SOC of carbon atoms and its large mobility^ make graphene a 



viable candidate for future spintronics applications, as demon- 
strated in non-local spin valve experiment a'^''^ . Spin trans- 
port studies yield Ti either by measuring the spin diffusion 
length As = l/\/d ■ vpy/TiT (where d = 2, 3 is the dimen- 
sion) or by a Hanle spin-precession experiment-. The first 
spin transport studies indicated a very short Ti w 100 ps spin- 
relaxation time which would be prohibitive for applications. 
We note that recent spin-transport experiments found a longer 
Ti w 2 — 6 ns for bilayer graphene which approaches the limit 
of applicabilit y'^i'^ . However, the experimental situation and 
the value of Ti is debated'^, as well as the appropriate theo- 
retical framework for the description of the experimentsi^^^. 

An alternative to measure the spin-relaxation time is con- 
duction electron spin resonance (ESR)-?. This yields Ti di- 
rectly from the homogeneous ESR linewidth, AB through 
Ti = I/7AB, where 7/27r = 28.0GHz/T is the electron 
gyromagnetic ratio. An advantage of the ESR method is its 
contactless nature, which also allows studying powder and air 
sensitive materials. It is however limited by the relatively 
large amount of samples required, which has yet hindered 
ESR studies on graphen©2^. 

It was recently shown by angle resolved photoemission 
studies that alkali intercalated graphite is a model system of 
biased graphene since the alkali atoms effectively decouple 
the graphene layersai^. Graphite intercalation compounds 
(GICs) have been known for decades^ and also, ESR was 
reported on them^S but to our knowledge no thorough study 
has been performed in order to unravel the relaxation mech- 
anism and in particular to study its anisotropy. Progress in 
the field of spin-relaxation in metals now allows a characteri- 
zation using the empiiical Beuneu-Monod plot"*'-, which tests 
the Elliott-relation (Eqs. ([T]i and (|2]i) based on empirical mea- 
surables alone. ESR studies on the GICs might also provide 
clues about the relaxation mechanism in biased graphene. 

Here, we report temperature dependent conduction electron 
spin resonance measurements on KCs compounds made of 
powder and HOPG (highly oriented pyrolitic graphite) sam- 



pies. KCs is the so-called stage 1 compound with the high- 
est available K intercalation level. We find that the momen- 
tum and spin relaxation times follow a similar temperature 
dependence, even if they are non-linear at high temperature. 
This proves that the Elliott- Yafet spin mechanism is valid in 
KCs and allows to place this material on the Beuneu-Monod 
plot. The measurement on the HOPG sample shows a size- 
able g-factor anisotropy, however it is mostly related to the 
anisotropy of the macroscopic susceptibiUty and the intrinsic 
g-factor is nearly isotropic. 



II. EXPERIMENTAL 

Stage 1 potassium doped graphite samples were prepared 
from round discs of HOPG (3 mm diameter, 50-70 /im thick- 
ness, mosaic angle; 0.4° ± 0.1°, "Grade SPI-1 HOPG", SPI 
Supplies Inc.) and graphite powder (grain size; 5-20 ijxn, 
Fisher Scientific Inc.). Prior to intercalation, graphite sam- 
ples were annealed at 500 °C in vacuum. Afterwards, sam- 
ples were handled in a Ar-filled glove box to avoid exposure 
to oxygen and water Doping was achieved through the two- 
zone vapor transport intercalation method^^ at 250 °C for 20 
hours with a temperature gradient of 5 °C. The color change 
of HOPG from gray to gold attests a successful doping to the 
KCs stoichiometry. No such significant color change is appar- 
ent for the powder samples probably due to surface roughness. 
The resulting samples were transferred to an ESR quartz tube 
and sealed under 20 mbar He for the measurements ^°. Powder 
samples were mixed with dilute Mn;MgO (1.5 ppm) to allow 
efficient microwave penetration. The ESR-silent and doping 
insensitive MgO separates the graphite grains, while Mn^+ 
has the added benefit of being a g-factor standard. Its g-factor 
is (7(Mn^"'")=2.0014^ and second order hyperfine interaction 
effects of the Mn^+ (Ref. [321) were taken into account to de- 
termine the g-factor of our samples. 

Experiments at 9 GHz (X band) were carried out on two 
commercial ESR spectrometers covering the 50-500 K tem- 
perature range with a typical microwave power of 10 mW and 
modulation of 0. 1-0.2 mT. 



III. RESULTS AND DISCUSSION 

A. The ESR llneshape 

The X-band spectra of the stage 1 potassium doped graphite 
in different sample forms are shown in Fig. [1] along with a 
photograph of the sample inside a quartz tube. Asymmetric 
derivative Loretzian lines, known as Dysonian lineshapes22iS 
are observed. The appearance of the Dysonian lineshapes 
is due to two effects. First, penetration of the exciting mi- 
crowaves, which is characterized by the penetration (or skin) 
depth (5) can be smaller than the sample size (d) for well 
conducting samples. Second, conduction electrons are mobile 
and diffuse while carrying along their spins even into regions 
of the sample which are not accessible for the microwave ex- 
citation. The diffusion is characterized by Tu, which is the 
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FIG. 1 . (Color online) X-band ESR spectra of KCg for different sam- 
ple forms. Note the asymmetric lineshapes characteristic for conduc- 
tive samples. Inset shows a gold colored KCg sample prepared from 
HOPG inside a quartz tube. 



time it takes for an electron to diffuse through the penetration 
depth. The Dysonian lineshape takes a number of different 
forms depending on the values of (5, d, Ti, and To. 

We found that the ESR data can be fitted with a mixture of 
Lorentzian absorption and dispersion curves with an almost 
1 ; 1 ratio^. This situation is the so-called "NMR limit" of the 
Dysonian lineshape and is described by Eq. (3.7) in the work 
of Feher and KipSi.. This occurs when To > Ti i.e. the elec- 
trons are diffusing through the penetration depth very slowly. 
This behavior in doped graphite was explained by Walmsley 
et al. for n-doped GICs [Ref. |35|l- The argument considered 
the conduction anisotropy of graphite, i.e. pc ^ pat and the 
platelet-like structure of HOPG such that sample dimensions 
in the ab crystalline plane are much larger than that along the 
c-axis. The conduction anisotropy and the good ab plane con- 
ductivity results in a relatively small penetration depth. How- 
ever, somewhat counterintuitively, the relevant diffusion term, 
which characterizes Tu is that along the c-axis. The diffusion 
being limited along this direction, the To ^ Ti holds which 
explains the experimental observation. 

The fact that we observe a similar lineshape for both the 
doped HOPG and graphite powder samples indicates that this 
situation also holds for the latter For all spectra, the asym- 
metry is moderate enough to fit the acquired spectra with such 
Lorentzian combination curves. The parameters of these fits 
yield the intensity, resonance field, and linewidth (HWHM; 
half width at half maximum of the Lorentzian resonance 



curves). 



B. The temperature dependent ESR linewidth 

The temperature dependence of the ESR spectra was mea- 
sured for the different sample forms. The extracted linewidths 
are shown in Fig. |2] Doped HOPG samples were studied 
only up to 450 K as for higher temperatures the linewidth 
changes irreversibly, possibly due to the evaporation of potas- 
sium atoms from the surface. The anisotropy of the linewidth 
was also studied, for an external magnetic field, Bq, parallel 
{Bq\\c) or perpendicular (Bq -L c) to the c-axis. 
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FIG. 2. Temperature dependence of the ESR linewidth in KCs (left 
axis) at 9 GHz. HOPG: Bo ± c (•) and Bo || c (o), powder (■<). Data 
on HOPG with Bo || c (if) from Ref. |23 is shown for comparison. 
Note the small linewidth anisotropy as found for the HOPG sam- 
ple. Dashed curve shows in-plane resistivity data (right axis) from 
Ref.lSa. with an axis offset in order to shows that it follows the same 
temperature dependence as the linewidth. 



The temperature dependence of the linewidth shows the 
same trend irrespective of the sample form. The curves 
mainly differ by a constant linewidth term: the powder sam- 
ple has a characteristically larger linewidth, which is either 
due to the larger effective suiface which gives rise to addi- 
tional relaxations or due to a higher concentration of impu- 
rities. In general, the linewidth can be written in analogy to 
the Matthiessen's rule for resistivity, i.e. as a sum of different 
contributions: 



AB(T) = ABhom(r) + ABi 



hom.O 



(3) 



Here, AShom.o is a temperature independent (i.e. residual) 
homogeneous relaxation such as caused by the surface (or im- 
purities) and ABhom{T) is the temperature dependent homo- 
geneous relaxation, which is in the focus herein. 

The linewidth data for the HOPG in the two orientations 
show a tiny anisotropy. The linewidth data for the HOPG are 



TABLE I. g-factor and in its shift with respect to go = 2.0023 for 
the different KCg compounds. The procedure to obtain the corrected 
g-factor is explained in the text. Results from Ref . |2^ are also shown. 
Errors are estimated from the variance of the data for different sam- 
ples. 
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2.0013 
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HOPG (measured) 


Bo-Lc 


2.0028 
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Bo II c 
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HOPG (Ref. 28.) 


Bo II c 
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HOPG (Ref. 28.) 
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7 



in good agreement with that reported by Lauginie et al. in 
Ref. |2^ which means that not only the temperature depen- 
dent AB\iom{T) is the same but there is a generic ABhom.o 
linewidth which is the same for HOPG samples from different 
sources. 

Our temperature resolution allows to observe a non- 
linearity in the temperature dependence of the linewidth 
which has been not reported, yet. It is best shown by a direct 
comparison to the temperature dependent in-plane resistivity, 
p{T) as it is done in Fig. [2] The two kinds of data are scaled 
together and are shown with an offset since the residual and 
temperature dependent contributions to AB{T) and p{T) do 
not necessarily scale together 

The temperature dependence of the in-plane resistivity, 
p{T), is quadratio^^i^ and is shown in Fig. |2] from Ref. |3g. 
For usual metals, p{T) is linear for temperatures T > 0d/2. 
The Debye temperature is 8d = 235 K in KCg^, whose value 
underlines the unusual nonlinearity of p{T) in KCg. It was ar- 
gued that this surprising quadratic temperature dependence of 
p{T) in KCs arises from electron-electron interaction^. We 
demonstrate in Fig. |2], that the temperature dependence of 
the ESR linewidth follows closely that of the resistivity. This 
leads us to suggest that the Elliott- Yafet theory is the appro- 
priate description of the spin-relaxation in KCg- 



C. The g-factor In KCg 

In the following, we discuss the value of the g-factor in 
KCs, which allows to follow the procedure of Beuneu and 
Monodi in verifying the validity of the Elliott-relation in KCs 
quantitatively (Eqs. ([T]i and (O). Room temperature 5-factors 
were obtained by comparison to the Mn^+ reference and the 
results are summarized in Table U along with data from pre- 
vious studies-**. While the present data for the HOPG show 
the same anisotropy as in the literature, an overall differ- 
ence of 2 — 3 • lO^** is observed for both orientations, whose 
origin is unclear™. We checked the consistency of the g- 
factor values for the HOPG and the powder samples with a 
numerical comparison: the anisotropic HOPG g-factor data 



was used to simulate powder spectra with a uniaxial g-factor 
anisotropy, which was then fitted with a Lorentzian. It yielded 
Ag = 1 • 10^^, which agrees with the experimental data for 
the powder as shown in Table J] 

A correction is required to determine the intrinsic g-factor 
in KCs due to its macroscopic magnetism. Again, usual al- 
kali metals (which are the textbook example of the EY the- 
ory) do not have a sizeable macroscopic magnetism. Graphite 
has unusual magnetic properties, as it has a relatively large 
and anisotropic diamagnetic susceptibility. This magnetism 
was associated with the orbital currents in graphitei22i and it 
not only affects the macroscopic properties but it also couples 
to the microscopic measurables. GICs also exhibit a sizeable 
macroscopic susceptibility but it is paramagnetic and its ab- 
solute value is smaller than that of pristine graphite. Such 
a magnetism is known in general to interact locally with the 
conduction electrons^, thus the measured 5-factor, ^meas, is 
different from the intrinsic one, gimr. This situation is well 
known in NMR, where local magnetism, e.g. those due to or- 
bital currents couple to the nuclei differently than the macro- 
scopic magnetismii. 

The resulting 5-factor in KCs is determined by the static 
susceptibility, xo, and a mean-field like coupling constant A 
and reads: 



5inu- (1 + ^X0), 



(4) 



2.0028 2.0030 




2.0010 



2.0015 



2.0020 2.0025 

g-factor 



2.0030 



2.0035 



FIG. 3. Illustration of g-factors in KCs: the go = 2.0023 free elec- 
tron p-factor (O)' data from Ref. 12^ (<1 and O), current measure- 
ments ( and ), and current measurements corrected by coupling 
to macroscopic susceptibility (M and ►). 



We require some simplifying assumptions to determine the 
intrinsic 5-factors in KCg: the coupling constant A is isotropic 
and is not affected by the alkali doping. We obtain A = 
-1.12 • 103 ^ from xr"''(^o II c) = -21.1 • IQ-^ sm 

(Ref.© and A.gS'^^Ph"'=(Bo || c) = +0.0473 (Ref. El. Here, 
we assumed that the whole g-factor shift in pristine graphite 
stems from the diamagnetic susceptibility and that its intrinsic 
g-factor is go. 

KCs has macroscopic susceptibilities of Xo'^^(-Bo II c) = 
1.02 • 10-6 £|H and Xo'^^B^ -L c) = 0.28 • 10"^ Sffiii 
(Ref. l42h . When these values are corrected by the diamag- 
netic core shielding and the Pauli term stemming from con- 



duction electrons^l, the true contributions from orbital mag- 
netism is obtained: Xo±^iBo \\ c) = 0.769 • 10"^ sm and 

Xo^\Bo -L c) = 0.029 • 10-6 ^. These values lead to the 
corrected, intrinsic g-factors in Table H] which are both posi- 
tive in sign and have a much lower anisotropy than the non- 
corrected values. The small g-factor anisotropy is in agree- 
ment with the similarly small anisotropy of the ESR linewidth. 



D. The microscopic theory of spin-relaxation in KCg 

The magnitude of anisotropy is important for the micro- 
scopic theory of spin-relaxation in intercalated graphite. This 
is also thought to be relevant for graphene as the GICs were 
suggested to act as its model system2ii2^; namely it was shown 
that the conduction band in KCs consists of Dirac cones (such 
as in graphene) with a Fermi surface separated by 1.35 eV 
from the Dirac point25,. Therefore, we attempt to adopt the 
description which was used in biased graphene to explain the 
spin-relaxation behavior in KCs. 

The kinetic energy in graphene is given by the low-energy 
Dirac Hamiltonian, which is obtained by expanding the tight- 
binding Hamiltonian around the corners of the hexagonal Bril- 
louin zone (the K and K' points), asii 



^Diiac = Vf {OxPx + CTyTzPy) - M, 



(5) 



where cr's describe the two sublattices of the graphene honey- 
comb lattice, vp = 10^ m/s is the Fermi velocity, r^ = 1 (-1) 
describes the K (K') degeneracy, and the /^j chemical potential 
accounts for a finite doping. For pristine, undoped graphene, 
fj. = 0. 

The resulting energy spectrum reads as a function of the 
momentum p as e(p) = wf|p|- This has to be supplemented 
with a spin-orbit coupling term to explain the spin-relaxation. 
It was argued previously that the intrinsic, i.e. atomic spin- 
orbit coupling is responsible for the spin-relaxation and g- 
factor shift for both graphene and the GIC&2J. provided the 
linear energy dispersion applies. The spin-orbit Hamiltonian 



is: 



i?so = A 



intr ' z'-' Z'^z: 



(6) 



where Ainn is the strength of the intrinsic SOC, Sz accounts 
for the physical spin, and cr^ measures the energy imbalance 
between the two sublattices. The above form of H50 satisfies 
time-reversal invariance, as expected. 

Eqs. ^ and ^ predict a significant anisotropy: the spin- 
relaxation would be sizeable when B L c and it would vanish 
for -B II c, since only s^ is involved in the above Hamiltoni- 
ans. On the other hand, the g-factor shift would be finite when 
B \\ c and zero when i? _L c. As mentioned, the experimental 
data does not show a sizeable anisotropy for either of these p a- 
rameters. At present, we do not have a consistent explanation 
for the difference between the theoretical model and the ex- 
perimental data. A possibility is that the simplest linear band 
model is not sufficient and that additional bands and spin-orbit 
couplings are involved such as e.g. those of the K+. Although 
these bands are well separated from the conduction bands!, 



they might influence the spin-relaxation properties due to a 
weak hybridization of the alkaU and graphite bands. A similar 
effect was invoked to explain the alkali dependent linewidth 
in alkali doped fullerides^. 

We note, that the lack of such extreme anisotropy was also 
observed in graphene spin transport measurements^. Interest- 
ingly, even the magnitude and the sign of the present ~ 25 % 
linewidth anisotropy agrees with the ~ 20 % value observed 
by Tombros et al. in Ref.l45l 



The spin-relaxation mechanism in KCg and the 
Beuneu-Monod plot 
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FIG. 4. The Beuneu-Monod plot showing the connection between 
(Ag)^ and ratio of the linewidth and resistivity (A_B/p) for pure 
metals (•) from RefQand for KCg for ± (■<) and || (►) orientations 
(error bars are smaller than the symbols), a) The original version 
with cgs units, b) the corrected plot with the dimensionless reduced 
linewidth {^i^B j psoujl-^ and considering the variation of Wpi among 
the materials. Solid lines correspond to Q.i/a2 = 1 and the dashed 
curve is the best fitting ai/ai = 10 as found in Ref.lj. 



in KCs. The Elliott-relation can be rewritten in terms of these 
experimental measurables as: 



AB 

P 



7 Ti 



7 (4 



i^gf 



(7) 



where Eq is the dielectric constant, cupi is the plasma frequency. 
In the summarizing work of Beuneu and Monod"* for most 



metals, -^^ ratio was found to be linearly proportional to 
(Ag)^ with a constant coefficient of 10^^ G/51cm (when us- 
ing cgs units), neglecting the variation of Wpj from metal to 
metal. 

The original Beuneu-Monod plot is shown in Fig.|4^. along 
with the present data for KCs, calculated from the present 
ABhom(r) and Ag data and from the p{T) in Ref.lM Cleariy, 
the data points for both HOPG orientations lie at an order of 
magnitude lower value than the majority of the other mate- 
rials. This can be explained by the relatively low value of 
ujp\ ~ 2.35 eV in KCs^ as compared to the usual value of 
> 5eV in elemental metals. In the original work of Be- 
uneu and Monod, variation of cjpi among the metals was not 
taken into account. We note that Petit et al. also suggested 
previouslyi^ that the lower utpi explains why data for alkali 
doped fullerides also fall relatively low on the Beuneu-Monod 
plot. 

In Fig. I^^., we show the corrected data by introducing the 
dimensionless reduced linewidth, jAB/peoUJpi, as a function 
of (Ag)^. With this correction, the data for KCs lie on the 
ai/a2 = 10 straight line, which was found to best fit the ex- 
perimental data for most elemental metals in Ref. Q We note 
that with this correction, the data for Rb and Cs does not seem 
to agree with this straight line in contrast to the original asser- 
tion. Whether Rb and Cs represents an anomalous situation 
herein or the somewhat old experimental results need to be 
revisited, requires additional work. Nevertheless, the present 
result confirms that KCs also follows the Elliott- Yafet theory 
of spin-relaxation and even the ai 2 parameters, which are 
sensitive to the band structure, are similar to that in ordinary 
elemental metals. 



IV. CONCLUSIONS 

In summary, we found the Elliott- Yafet theory of spin- 
relaxation to be valid for the KCs stage 1 graphite intercala- 
tion compound, with the proportionality of the homogeneous 
linewidth and in-plane resistivity in agreement with the g- 
factor shifts as for most elemental metals. It remains, how- 
ever open for further investigations whether the result can be 
applied directly for the spin-relaxation in biased graphene as 
expected based on the band structure calculations and ARPES 
studies2ii2^. 



The above discussed temperature dependent ABhomiT), 
the p{T) data from the literature (Ref. 1361) (both are free from 
the temperature independent terms) and the corrected value 
for Ag allows us to discuss the validity of the Elliott-relation 
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